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A simple analytical model has been developed for the electric and ponderomo 
tive (trapping) potentials in linear ion traps. This model was used to calculate the 
required voltage drive to a mercury trap, and the result compares well with exper- 
iments. The model gives a detailed picture of the geometric shape of the trapping 
potential and allows an accurate calculation of the well depth. The simplicity of 
the model allowed an investigation of related, more exotic trap designs which may 
have advantages in light-collection efficiency. 


I. Introduction 

Radio frequency (RF) quadruple ion traps have great 
importance to the development of new atomic frequency 
standards and high-precision measurements. The three- 
dimensional quadrupole trap has mainly been used. The 
ideal geometry for the electrodes of this trap is hyper- 
boloids of revolution, which produce a pure quadrupole 
field. Practical traps are built with a different, more open 
geometry because of the need to collect light efficiently 
from the ions. Ion fluorescence is used to determine its 
quantum state, and the collection efficiency is an important 
factor in the signal-to-noise ratio. Deviations in the elec- 
tric field and trapping potential from a pure quadrupole 
are usually not considered in detail due to the difficulty in 
calculating them. 

In order to increase the number of trapped ions with- 
out degrading the frequency stability, the Time and Fre- 
quency Systems Research Group recently introduced into 
frequency-standard research a linear trap based on the 
quadrupole mass spectrometer [1]. A side benefit of this 


trap geometry is the ability to calculate the electric and 
trapping potentials to good accuracy with a simple an- 
alytic model. In this article, the model is developed and 
used to predict an ion resonance frequency, which has been 
experimentally measured. The article then shows how the 
model can also be used to investigate similar but more 
exotic geometries which may be advantageous in some ap- 
plications. 

The linear trap consists of four parallel cylindrical 
rods arranged with their centers on the corners of a square. 
An RF voltage is applied to the rods so that nearest neigh- 
bors have opposite polarity. This creates an alternating 
two-dimensional quadrupole electric field between the rods. 
The field confines ions along the center axis of the trap 
by the ponderomotive force, just as in a quadrupole mass 
spectrometer. Two endcap electrodes with a DC bias volt- 
age applied confine the ions axially. This article considers 
only the center axial region, far enough from the endcaps 
so that the field is essentially two-dimensional. The DC 
fields from the endcaps decay exponentially along the axis, 
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so that this restriction applies to the vast majority of the 
trapping volume. 

The electric fields due to four cylindrical rods have 
been studied for years by designers of mass spectrometers 
and particle accelerators. In order to achieve the closest 
approximation to a true quadrupole field, it has been found 
that the ratio of the rod radius to the distance between the 
rod center and the trap axis reaches its optimum value at 
0.5342, This value was first determined by measurements 
on a quadrupole accelerator magnet [2] and later repro- 
duced with numerical calculations [3]. 

A radius-spacing ratio of 0.5342 produces a geometry 
that is too closed to allow efficient light collection from 
ions. A pure quadrupole field is not really necessary if ion 
confinement is the only goal; hence, in this case the ra- 
tio can be reduced. The trap discussed here has a ratio 
of 0.25, and even smaller ratios might be desirable. The 
smallness of the ratio leads directly to a simple approxima- 
tion, namely, the fields produced by infinitely thin rods. A 
conducting rod with an applied voltage has some induced 
charge, and as the rod diameter is reduced, all the charge 
coalesces into a line. This idea can be used to transform a 
two-dimensional boundary- value problem into a much sim- 
pler calculation of the potential produced by fixed sources. 
All the calculations that follow are based on fields gener- 
ated by an array of uniform parallel line charges. 


II. Model for a Four-Rod Trap 


Changing to polar coordinates, it is found that 


|F | 2 = 


16r 2 

r 8 — 2 r 4 cos 4<£ + 1 


r GM) 


(i) 


where V = 16. One can expand G about the origin to 
obtain 


G(r, $) w r 2 1 + 2 r 4 cos 4<j> + r 8 (l + 2 cos 8<£) + • • ■] 


where the leading term r 2 is an isotropic harmonic po- 
tential, and the higher-order terms have at least fourfold 
symmetry. It turns out that any two-dimensional field con- 
figuration that vanishes at a point produces an isotropic 
harmonic trapping potential in lowest order about that 
point. In the trap configurations analyzed later, the func- 
tion G is always defined to have a leading term r 2 , and T 
is used for the numerical factor. 

The trapping potential model function G(x i y) is 
shown in a three-dimensional plot in Fig. 1, and shows a 
center well rising up to singularities at the rods. Ilalfway 
between the rods are saddle points beyond which the po- 
tential drops again. The height of the saddle points sets the 
maximum energy an ion can have and still stay trapped, 
i.e., the well depth. The well depth can be calculated by 
setting <f> = tt/ 4 in Eq. (1) and finding the maximum value 
of G(r, 7t/4). The maximum occurs at r s — 3~ 1 / 4 =: 0.760, 
and has a value of G s — 9/16\/3 — 0.3248, where the 
subscript s is the saddle point. 


Since the model field is two-dimensional, the method 
of complex variables can be used. This is not necessary, 
but makes the calculation of the trapping potentials a little 
easier. Adopting the notation of Landau and Lifshitz [4], a 
positive line charge at z = zq produces a complex potential 
w — — log(z — zo), whose real part is the ordinary electric 
potential. Dimensionless quantities are now used to cal- 
culate the geometric form of the potentials. Scale factors 
are introduced later to calculate real trap parameters. The 
model for the quadrupole trap consists of two negative line 
charges at z = ±i and two positive line charges at z = ±1, 
This produces the complex potential 



Now that the electric and trapping potentials pro- 
duced by four line charges have been calculated, a real trap 
may be modeled. The electric equipotentials of the model 
are plotted as contours in Fig. 2. Close to the charges, 
they have a nearly circular shape. The equipotentials of 
a real trap are exactly circular at the electrode surfaces. 
The deviation of the modeFs equipotential from circularity 
is calculated by finding its horizontal and vertical “diam- 
eters” and taking the difference. The electric potential V 
is the real part of and is given by 

7 v _ s 4 + y 4 + 2 (* 2 - s/ 2 + * 2 y 2 ) + 1 _ B 2 , 2 \ 

e x 4 + y 4 + 2(y 2 - x 2 + x 2 y 2 ) + 1 


The ponderomotive trapping potential is proportional to 
the square of the electric field F [5]. To calculate this, note 
that — F x + iF y — dw/dz ) so that 


The equipotential surrounding the charge at z = 1 is to be 
examined here, so first the intersection of the equipotential 
with the z-axis is found by using 


|F| 2 = 


dw 2 
dz 


16|z| 2 


z 4 - 1| 2 


11 + z 2 | 

11 - a; 2 | 


= D 


( 3 ) 
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For a fixed value of V, Eq. (3) has two positive and two 
negative solutions. If the smallest positive solution is de- 
noted xo, and the other positive solution is denoted xi, 
then Eq. (3) produces the following relationships: 

i [bTT 


2 

1 Xo ~ x 0 ~ Vfl 5 - 1 

/. = x ° + fi _ 1 + Xo _ B 

2 x 0 y/B*~- 1 

which are valid for B > 1. Notice that there is a symme- 
try between B and 1 /B corresponding to the symmetry 
between V and —V. The case of B < 1 applies to the 
region near the charges at z — ±*. 

Point xo corresponds to the inner surface of a rod, x\ 
to the outer surface of a rod, d to the horizontal diameter of 
the rod, and c to the center of the rod. The “rod” defined 
here is a hypothetical electrode shaped so that four of them 
will produce a field distribution equivalent to the four line 
charges. To calculate the vertical diameter, one can use 
the intercepts ±t/o on the line x = c * U s i n S Eq. (2) results 
in 


(4a) 

(4b) 

(4c) 


2 B (\/(5 2 + 2) - B) - 1 ! 

y 0~ B 2_x ~ 02 

The approximation is very good for B > 2, and improves as 
B increases. The fractional difference of the rod diameters, 
(d~ 2 yo)/d, scales as l/(2£ 2 ), and is consistent with the 
idea that the equipotentials become more circular as the 
line charge is approached. A curious feature of the model is 
evident in Eq. (4c). The center of the rod c depends on B 
and does not coincide with the line charge at x = 1. Seen 
another way, the line charge sits at the geometric mean of 
the inner and outer surfaces of the rod (xo and xi) and 
only approaches the arithmetic mean as B becomes large. 

So far, it has been assumed that the boundary is in- 
finitely far away. This is allowed since the potential V in 
Eq. (2) tends toward zero for large r. The fact that the to- 
tal charge is zero assures this. Denison [3] has numerically 
calculated the effect of a circular grounded boundary at 
r = 1.65 for an optimized quadrupole with a radius-spacing 
ratio of 0.5342. He found that the boundary changed the 
optimum ratio by about 1 percent. Although a more open 
trap may be more sensitive to a boundary, practical traps 


will have their boundaries at greater distances than 1.65. 
Therefore, boundaries will be ignored in this treatment. 

In order to use the model to calculate the parameters 
of a trap, the “squashed” rods of the model are associated 
with the round electrodes of the trap, and the linear di- 
mensions are scaled appropriately. Dimensions of the real 
trap are given in capital letters; model dimensions are in 
lower case. The distance of the line charges from the trap 
axis is A, and they will not be exactly at the center of the 
electrodes. The inner and outer x-axis intercepts of the 
trap rods are called A" 0 and X Xy respectively. Therefore, 
the dimensionless coordinates are given by 


II 

O 

N 

(5a) 

H 

II 

(5b) 

A 2 = X 0 Xi 

(5c) 


The mercury ion linear trap has dimensions Xo = 7. 62 mm 
and X x - 12.7 mm, implying A = 9.84 mm (compared to 
10.2 mm for the true center), xo = 0.774, x\ = 1.29, and 
B — 3.99. The horizontal and vertical diameters of the 
equivalent model trap rod differ by 3 percent, providing a 
good approximation to a circle. 

In order to test the predictive power of the model, the 
applied voltage needed for the mercury ion trap is calcul- 
cated to produce a natural ion resonance frequency u> of 
27r x 48.5 kHz. For small displacements, the leading r 2 
term in the trapping potential leads to harmonic motion. 
The ions are detected by amplifying the current they in- 
duce on the trap rods at the natural frequency, and the 
amplifiers are tuned to 48.5 kHz. The driving voltage is 
applied at a frequency Q of x 500 kHz in a balanced 
mode, so that the peak voltage on two opposing rods is Vo 
and the voltage on the other two rods is - Vo with respect 
to the vacuum system. The balanced drive keeps the trap 
axis at zero potential, which allows the application of only 
DC bias to the endcaps. The amplitude of the applied RF 
voltage is ramped until a resonance signal appears, indi- 
cating that the natural frequency of the trap matches the 
frequency to which the amplifiers are tuned. 

Calculating the required applied voltage uses the fact 
that on the inner surface of the rod Xo, the electric poten- 
tial V equals the applied voltage Vo, so that 

V _ In |(1 + g 2 )/(l — s 2 )l In |(1 + g 2 )/(l — g 2 )l 

Vo “ ln|(l + *2)/(l-*2)| lnS 
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and the electric field is 


and have the value G s = 1/4. The electric potential can 
be written as 


|£| 2 = 


rK 0 2 G(x,y) 

(A\nB) 2 


e 2V = x 4 + y 4 + 2(x 2 — y 2 + x 2 y 2 ) + 1 = B 2 (9) 


The ponderomotive trapping potential 'P is given by 


_ e 2 E 2 _ Te 2 VoG(x 1 y) 
4 mQ 2 4mQ 2 (A\nB) 2 


( 6 ) 


One can calculate the intercepts of an equipotential on the 
y-axis around the charge at z = i and obtain relations 
analogous to those of the four-rod trap: 


where m is the mass of the ion [5J. Using the harmonic 
approximation for G leads to an expression for the natural 

J/O = V 1 — B 

(10a) 

resonance frequency 

yi = V 1 + B 

(10b) 


VTeVo 

s/2 m£l 2 A 2 In B 


( 7 ) 


For mercury isotope 199 and the trap parameters given 
above, Eq. (7) predicts a peak drive voltage of 93.7 V. The 
experimental value is 100 ± 5 V, giving agreement almost 
within experimental error. The well depth can also be 
calculated using Eqs. (6) and (7) as 

*. = 2 =^! ( 8 ) 

yielding a value of 3.00 eV for the parameters here. 


III. Model for a Two-Rod Trap 

The simplicity of this model makes it useful for ana- 
lyzing other two-dimensional trap geometries. Any system 
that produces a point with a vanishing electric field can 
potentially trap ions. The simplest system consists of two 
rods driven by the same voltage with respect to a dis- 
tant boundary. This configuration has a zero field point 
midway between the two rods. A trap of this type was 
demonstrated with oil droplets by Straubel [6] very early 
on, and may be of use in frequency standards due to its 
wider viewing angle. The model for this trap is two equal 
negative line charges at z = ii. The complex potential is 

w = ln(z 2 + 1) 

and the geometric form of the ponderomotive potential is 

r 2 

r 4 + 2r 2 cos 2<f> + 1 


for B < 1. Combining Eqs. (10a) and (10b) to calculate 
the scaling law for a real trap results in 2 A 2 = Y 0 2 + Yj 2 , 
where Yo and Y\ are the inner and outer surfaces of the 
rod. For a trap with rods the same size and spacing as 
two opposing rods of the mercury ion trap, A = 10.5 mm. 
This time the line charge is more distant than the rod 
center (10.2 mm). A calculation of the diameter difference 
gives a deviation from roundness of 3 percent. 

In order to calculate the applied voltage, the surround- 
ing boundary must be taken into account, because the bars 
must be driven with respect to something. According to 
Eq. (9) and Fig. 3, the equipotentials for large r become 
roughly circular with a value of 2 In r, and it is assumed 
that the boundary follows one of these contours. The po- 
tential on the trap axis is zero and the potential on the 
rod is In B. The applied voltage will be proportional to 
2 In r — In B, where r = R/A and R is the radius of the 
surrounding boundary. The trapping potential can then 
be written as 


* = 


Te 2 V$G 

4mQ 2 A 2 (2 In r — In B ) 2 


The dimensions of the vacuum system correspond to a 
value of 2.5 for r. In order to duplicate the natural res- 
onance of 48.5 kHz with the two-rod trap, a peak drive 
voltage of 364 V is necessary, an increase by a factor of 3.9 
over the four-rod trap. The well depth is still governed by 
Eq. (8), and depends on the boundary only through the 
voltage necessary to maintain u). The well depth then is 
2.63 eV. The price paid for an increased viewing angle is 
a larger drive voltage and a slightly smaller well depth. 


with T = 4. Contours of the electric equipotentials and a 
plot of the ponderomotive potential are shown in Figs. 3 
and 4. The two saddle points are on the x-axis at r B = 1, 


An added complication with this trap is that the elec- 
tric potential at the trap axis is now oscillating with respect 
to the vacuum system at the drive frequency. This makes 
an added AC bias to the endcaps necessary to keep them 
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at a constant potential above the trap axis. Straubel did and for the outer rods (B 2 < 1) 
not use endcaps, but relied upon fringing fields from the 

ends of his rods for axial confinement. d 2 = B 


IV. Models for a Three-Rod Trap 

Although the viewing angle is larger in the two-rod 
trap, the trapping point is still between the two rods. It 
may be desirable for some applications to trap well out- 
side the electrode structure. This becomes possible with a 
three-rod trap. A simple example of a three-rod trap model 
consists of a positive line charge at the origin flanked by 
two equal negative line charges on the real axis at z = ±1. 
The resulting complex and ponderomotive potentials are: 



and 



Once again, the intercepts of the outer rod obey the rela- 
tions of Eq. (5). The boundary is again assumed to con- 
form to a long-range equipotential at dimensionless radius 
v with value In v . All the potentials can be offset by In r 
to keep the boundary at ground and obtain the following 
form for the applied voltages: 

Vi ln(B x /r) 

V 2 In (B 2 /r) 


\Ff = 


r 4 -f 2 r 2 cos 2<fi + l 
r 2( r 4 _ 2 r 2 cos 24> 4- 1) 


( 11 ) 


The electric field vanishes at z = ±i, and these are the 
trapping points. The electric and trapping potentials are 
shown in Figs. 5 and 6. 


If Eq. (11) is expanded about the trapping point z = i, 
then T — 4. The expression for G is too complicated to 
be very useful. The lowest saddle point can be found from 
Eq. (11) at r s = y/2 + y/E = 2.058, and <j> = tt/ 2 where 
G s = 0.0217. This gives a well depth of 200 meV for the 
parameters used here. The question of the applied voltage 
is slightly more complicated than in the previous two cases. 
Since the total charge is not zero, the boundary must again 
be included in the calculation. There are now three volt- 
ages involved: the central rod V\> the outer rods V 2l and 
the boundary V3. The ratios between these voltages can 
be calculated from the expression for the potential 


Using a value of 2.5 for r and the same rod diameters and 
spacing as before obtains peak applied voltages of 49.4 V 
and -214 V for V x and V 2 . A suitably tapped transformer 
could be used to provide these drive voltages. Alterna- 
tively, the ratio of the inner and outer rod diameters could 
be adjusted to force B\B 2 = r 2 so that V\ ~ —V 2 . Even 
with this simplification, an additional AC voltage is nec- 
essary to bias the endcaps as in the two-rod case. 

The influence of the boundary can be removed by 
choosing the total charge to be zero. This can be ac- 
complished by doubling the center charge. Unfortunately, 
there is no trapping point with this configuration while the 
charges are in a line. Moving the center charge down along 
the imaginary axis to z — —ib creates a trapping point at 
z = i/b. The case of b = 1 has been analyzed, and the re- 
sults are shown here. The electric and trapping potentials 
are plotted in Figs. 7 and 8. The scaling factor T is 1/4, 
and the saddle point is on the y-axis at y s = 1.839 with 
G $ = 0.0728. The electric potential is 


x A 4- y 4 + 2(y 2 - x 1 4- x 2 y 2 ) + 1 ^ q2 
x 2 4- y 2 


x A 4- y 4 4- 2(y 2 - x 2 4- a V) 4- 1 
(x 2 -f y 2 4- 2y -b l) 2 


and the various intercepts. For the center rod (B\ > 1), 
the relation between the B value and the dimensionless 
diameter is 



and the x intercepts of the outer rods still follow the re- 
lationships of Eqs. (4) and (5), with B replaced by l/B\ 
(for Bi < 1). The y intercepts of the lower rod obey the 
relationship below (for B 2 > 1): 

-B 2 ± y/2B 2 - 1 
2/0(1 B 2 - 1 
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and Eq. (5). Once again, the voltages applied to the rods 
will not be equal and opposite unless the diameters are 
adjusted so that B\ = 1/B 2 . The endcaps need only DC 
bias because the potential is zero at the trapping point. 

The small value of I\ which is 1/4, makes the drive 
voltage for the outer rods eight times that of the four-rod 
trap for the same ion resonance frequency. This fact, plus 
the small well depth, may limit the usefulness of this trap. 
In both three-rod traps, however, the low saddle point lies 
only on one side of the trapping region. The potential 
barrier on the opposite side is much higher. These traps 
may be useful for trapping macroscopic particles, where 
gravity plays an important role. The trap rods could be 


oriented horizontally so that the low saddle point lies above 
the trapping region. 

V. Conclusions 

This article has presented a simple model for linear 
ion traps that permits the accurate calculation of trapping 
parameters and gives a detailed picture of the potentials. 
The model was used to analyze some new trap geometries, 
and their advantages and disadvantages were discussed. 
Each trapping geometry is characterized by the two pa- 
rameters T and G,, which determine the dependence of its 
natural resonance frequency and well depth on the applied 
RF voltage and trap dimensions. 
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Fig. 5. The electric potential for a three «od I trap is shown .In this ^ elcctrIc polentia , for a , hr ee-rod trap with zero net 

contour plot. The range of potentials shown is -1 to 1.5 at Inter- ^8 jg 8hown , n , h|8 contour pIol ^ range of poten tlals 

va s 0 shown Is —2 to 4 at Intervals of 0.5. 
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Fig. 8. The ponderomotlve potential function |F*/4| for a three- 

Fig. 6. The ponderomotlve potential function | F 2| for a three-rod rod quadrupole trap with zero net I ’l"* ""o^ 

quadrupole trap is plotted in three dimensions. The singularities dimensions The singular! es are truncated at the value 0.1 to 

are truncated at the value 0.5 to emphasize the potential well at emphasize the potential well at z - /. 





